Is the supersolid superfluid? 

Dimitar I. Pushkaro^|*] 

Institute of Solid-State Physics, 
Bulgarian Academy of Sciences, 

Sofia 1784, Bulgaria and 
Fatih University, Istanbul 3400 

An analysis of previous theories of superfluidity of quantum solids is presented in 
relation to the nonclassical rotational moment of inertia (NCRM) found first in Kim 
and Chan experiments. A theory of supersolidity is proposed based on the presence 
of an additional conservation law. It is shown that the additional entropy or mass 
fluxes depend on the quasiparticle dispersion relation and vanish in the effective mass 
approximation. This implies that at low temperatures when the parabolic part of the 
dispersion relation predominates the supersolid properties should be less expressed. 

PACS numbers: 67.80.Mg 



Introduction 

The experiments of Kim and Chan |l| breathed new life into the old idea of possible 
superfluidity of solids. A quantum solid possessing superfluid properties has been called a 
supersolid. Originally a supersolid should be a crystalline body where a nondissipative mass 
current can occur. This should correspond to the superfluid state of liquid helium (helium-II) 
observed by Kapitsa, and explained theoretically first by Landau. The superfluidity is now 
well studied and a number of effects has been found, predicted and explained. Between them 
is the change of the rotational moment due to the fact that the superfluid fraction cannot 
be involved into rotation at velocities less than the critical one. The qualitative explanation 
of such a behavior according to Landau is that at small velocities no excitations can be 
generated. A successful hydrodynamical description is given by the so-called two-fluid (or 
two-velocity) hydrodynamics. From a mathematical point of view the new element in the 
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two- velocity hydrodynamic equations is the potentiality of the superfluid velocity which 
reads Vg — V/i with for the chemical potential in the frame where = 0. As a result, 
a new vibrational mode, the second sound, appears. The phase transition into a superfluid 
state is well defined and the corresponding changes of the thermodynamic characteristics 
are well investigated. 

The quantum-mechanical consideration connects the superfluidity with the Bose-Einstein 
condensation (BEC). Later on, such kind of condensation in the momentum space was 
observed in some gases as well. This is the reason to talk on a macroscopic quantum state 
described by the condensate wave function. 

In their works Kim and Chan have observed a nonclassical rotational moment (NCRM), 
i.e. a rotational moment of inertia which changes its value with temperature in a way bearing 
a resemble with the Kapitza experiments with rotating liquid helium. They argue that this 
is enough to conclude that the body has been in a supersolid state and that the superfluidity 
has been finally observed in all three states of matter (gas, liquid and solid). Lately, the 
term supersohd has become a synonym of a body with NCRM. 

The first reasonable question is whether the NCRM implies superfiuidity (supersolidity) . 
Is the supersolid state "superfiuid" or this is an evidence of a new phenomena, maybe more 
interesting and famous than superfluidity, but nevertheless, of different kind. The existing 
experimental observations and theoretical analysis did not give an unambiguous answer yet. 

Originally, the concept of supersolid appeared for a crystalline body inside which a nondis- 
sipative (macroscopic) mass current can exist. First considerations (Penrose and Onsager, 
Andreev and Lifshits, Legget, Chester etc.) had a crystalline bodies in mind. Defects in 
such crystals are imperfections of the crystal lattice, or lattice with an ideal periodicity but 
with a number of atoms less than the number of the lattice sites ( Andreev-Lifshits) . The 
flrst question is therefore if the experiments can be understood from such a point of view. 
Most probably this is not the case. 

Let us first consider the validity of the Landau derivation of the critical velocity. In liquid 
helium, the energy in the frame where the superfiuid velocity is zero can be written in the 
form: 

E = Eo + Pqv + ^Mv^, P = Po + Mv (1) 



3 



The same relation for an elementary excitation e{p) reads 

E = e{p) + pv + ^Mv'^ (2) 

where p is the momentum in the frame where v = 0. The least possible change in energy 
due to the excitation created is e{p) — pv and should be negative in order to reduce the 
energy of the system. This yields 

V > e{p)/p. (3) 

It is worth noting that equation ([T]) is always valid because it follows directly from the 
Galilean principle for macroscopic quantities. Relation ([2]) corresponds to the microscopic 
characteristics of the elementary excitation. In a homogeneous and isotropic (Newtonian) 
space these two relations coincide. However this is not the case in a crystalline solid where 
quasiparticle states are classified with respect the quasimomentum, not momentum. Quasi- 
momentum is simply a quantum number which apears due to the periodicity of the lattice. 
There are no Galilean transformations for quasiparticle characteristics. The transformation 
relations which replace the Galilean ones were derived in The macroscopic 

momentum (the mass flux) is not a mean value of the quasimomentum, but of the product 

de 

m— with k for the quasimomentum. In addition, phonons in crystals have zero momentum 

OK 

and do not transfer mass in contrast to the phonons in liquids. All this implies that the 
Landau criterion ([3]) does not work in crystalline bodies. Its 'quasiparticle' analogue should 
look like 

, mde dine 

or 

dine , , 

mv < (5) 

But, there is still a question what is, say, the phonon qusimomentum in the co-moving 
(with the superfiuid fraction) frame. In addition, m = for acoustic phonons. To avoid any 
misunderstanding, let us stress again that whatever the dispersion relation of the elementary 
excitations and the spectrum classification parameter (momentum or quasimomentum), the 
macroscopic fluxes have to obey Galilean relation ([1]). 

Next, it is very important that the conservation laws (which are the basis of the hydrody- 
namics) exist only in an inertial laboratory frame. And this laboratory frame is privileged, 
not Galilean (see more details in Ref. jsl). 
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If one considers Bose-condensation of quasipartides then such a condensate in a crystal 
must be characterized by a value of the quasimomentum in a privileged frame. 

Finally, the particles or quasipartides (say vacancions) undergoing Bose-condensation 
should interact weakly enough. It was shown [6] that the vacancion gas is most ideal near 
the middle of the energy band, not in the bottom. 

It is seen, therefore, that the situation in a crystalline body is completely different com- 
pared to liquids and gases. 



Nevertheless, the first hydrodynamical theory of the superfluidity of solids |8| was devel- 
oped in a close analogy with the Landau theory of helium II. Andreev and Lifshits introduced 
two velocities for a normal and a superfiuid fraction of the solid and apply the potentiality 
condition for the superfiuid velocity, = V/x. They used the Galilean invariance, so in the 
frame, where the superfiuid component is in rest (v^ = 0), the energy per unit volume is: 

E = pVs'^/2 + p.Vs + £, j = pv^ + p (6) 

where j is the momentum per unit volume equal to the mass flow while p is the momentum 
in the frame with = 0, 

e = e{S,p,Wik) (7) 

is the internal energy as a function of the entropy, density and the distortion (not symmetric) 
tensor 



dui 



dxk 

The tensor of small deformations is as usually equal to 

I ( dui ^ duk 
Uik = - < -7{ h 



2 [ dxk dxi 

and its trace equals the relative variation of the volume 

Uii = Wii = 6V/V 

A new point is that this trace is now not connected to the density variation with the known 
relation, i.e. 

wu ^ (8) 



In this notation 



de = TdS + XikWik + i^dp + (v„ - Ws)dp. (9) 



A standard procedure follows based on the conservation laws: 

S + div (5v„ + q/T) = R/T, {R > 0) (11) 

V, + V<^ = 0. (12) 

The unknown quantities Ilik,(p,(l, R have to be determined so as to satisfy the redundant 
energy conservation law: 

^ + divQ = 0. (13) 

The time derivative of E reads: 

E = TS + K,^^- /idiv j - div (^j) + j (14) 
/• \^ 9Ilik . 

- [i- P^n)y'^-Vni^ hV^VgdlVJ 

OXk 

= -div + ST\ri + v„(v„p)^ + T{S + divS'v„) 

+ -^ik-g^ + (j - PV„) V I ~ ~^ ) ~ PV„V/i 

d 

- Vni-^ {Ilifc - pVsiVsk + VsiPk + VskPi 

OXk 

+ [-£ + TS'+ (v„ - Vs)p + yup](5ifc} -yudivj. 

Here, a term of the form VniXki^r~' neglected as cubic in "normal motion". With the 

dXi - 

aid of conservation laws the time derivative of energy was written in the form 
E + div <(^- + + STvn + v„(v„p) + q + (f{i- pv„) + 



This yields the following expressions for the fluxes: 



n^A; = pVsiVsk + VsiVk + VnkVi 

+ [-£ + TS'+ (V„ - Vs)p + - Aifc + TTifc, (16) 



if = ^ + p + 7/>. (17) 
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Q 



+ /i j + ST\n + v„(vnp) + q 



and the dissipation function of the crystal is 



R 



V^div (j - pwr. 



qVT 



[Vnk - Uk 



ik 



(18) 



(19) 



dxu ^-^ T ' ' dx- 

We shall not write here the relations between iTikjipyq and (vn — u) that follow from the 
Onsager principle and the positivity of the dissipative function. The main consequence is 
that with neglecting dissipation one has v„ = u. The normal motion is therefore the motion 
of lattice sites (which may not coincide with given atoms). The superfluid flow could, hence, 
be possible at a given (even not moving) lattice structure. 

However, instead of (1151) the time derivative E can also be written in the form: 

2 



E + div 



+ yu j + STvn + v„(v„p) + q + V'(j - pv„) + 



+ VnkT^ki — -KkVnk 

= R + Tlik— h V'dlV (j - pV„) + -—- + XikT^[Uk - Vnk), 

OXi I OXi 



which leads to other expressions for fluxes. 
In this case the nondissipative theory yields: 

2 



E + div 



+ Vnk 



^ + /i j j + STVn + V„(v„p) + ((^-/i-y )(j - pV„) + 



Ilfci - pVsiVsk + VskPi + VniPk- 

[-e + TS + (v„ - Vs)p + np]Sik 

I - pVsiVsk + VskPi + VniPk- 



[-e + TS + (v„ - Vs)p + pp]6, 



ik ' 



dv, 



nk 



dxi 



+ 



duk 

+ {ip- p - y )div (j - pv„) + \k-Q^ 



and hence. 



^ik = PVsiVsk + VsiPk + VnkPi + 
,2 



-e + TS + (v„ - Vs)p + pp- 



Hk 



2 



A 



ik 







(20) 



(21) 



(22) 
(23) 
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The procedure used is, therefore, not unique. The relation v„ = u was not derived, but 
presupposed. In fact, the consideration started as a three-velocity theory (u,Vn,Vs) and 
the identity li = Vn follows from a condition the time derivative of the total energy be not 
dependent on the which is not well grounded. Next, the conservation laws are written 
in the system where Vg = and this is not the laboratory frame in which the lattice cites 
are in their equilibrium positions. 

That is why we turned to another approach based on our theory of the quasiparticle 



kinetics and dynamics in deformable crystalline bodies ^, 0, Ej. This theory works with 
exact (in the frame of the quasiparticle approach) selfconsistent set of equations including 
the nonlinear elasticity theory equation and a transport Boltzman-like equation valid in the 
whole Brillouin zone of quasiparticles with arbitrary dispersion law. The theory is developed 
for crystalline bodies subject to time-varying deformations and arbitrary velocities. 

Partition Function and Thermodynamic Relations 

Let us consider a gas of quasiparticles with dispersion law £:(k) at low temperatures, 
when the frequency of normal processes is much larger that of the Umklapp processes, i.e. 

The distribution function nfc(k, r, t) corresponds to Smax with conserved energy E, quasi- 
particle density n, quasimomentum K and momentum (mass flow) j defined, respectively, 
as: 



Sir,t)= I sHdk, (24) 
where s[nk] = (1 + n^) ln(l + n^) — Inn^ 



E{r, t) = J EkUkdh (25) 

n(r, i) = j rikdk (26) 

K(r, t) = jkukdk (27) 

Uk dk, dk = dkidk2dk3 (28) 



This yields 



. - V.k-mW.((9£fe/(9k) -/i . 
nfe(k, rt) = <i exp ( ) - 1 



with V, W and /i for Lagrangian multiphers. Varying S yields 

T6S = SE- Y.6K - W.^j - ^i6n 
n^E-TS- V.K - W.j - 

and respectively 

dil = -SdT - K.dV - ydW -ndn 
The nondissipative equations involve the following conservation laws: 

n + div J = 0, J—j/m 
dji. ^ dUik ^ Q 

dt dxk 

dKi dLju 

— 1 ^ — m = 0, 

dt dxk 
S + divF = 0, 

E + div Q = 

To second order with respect to velocities V and W one has: 

Ji = nVj + UijWj, 

where 



dkidki 

The local drift velocity U = j/mn° is then 



-1 



The mass flux is, therefore, not coUinear to any of velocities V and W. Analogously, 

Ki ^ PiiVi + mn'^Wi, 
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To second order in velocities the diagonal terms ofthe quasimomentum flux tensor coincide 
with the thermodynamic potential fi(T, V,W,/x) |3|, y, [Sj: 

Lij = J ki^-Uk dk = 6ij (38) 

and the momentum flux tensor has the form: 

= -nu = T j ln(l + n°) z/i;(k) dk (39) 

The energy flux is: 

Q,{y,t) = j ek^nkdk = W'V, + {TS,i + fin,i)Wi (40) 

where 

Su = j s[nl]u,idk (41) 

and = E° - is the enthalpy at V = W = 0. Hence, W^Y is the energy flux known 
from the classical hydrodynamics, and there are additional terms due to the supersolid 
behavior. The full hydrodynamic system consists then of four equations: 

n + ndivU = 0, (42) 

mnif, - ^ = (43) 

OXi 

P^s^ - + n\6u - f3u)m = (44) 

E + iy°div V + Ts(^-'^)^ = (45) 



S n J dx. 
where 

I3ii = P^knkl/{n^f (46) 
Taking into account the thermodynamic identity: 

dE = TdS + pdn + \.d'K + W.d] (47) 

the energy conservation law can be replaced by entropy equation: 

^ + ^divU + ^f%-^')^ = (48) 
\ D n J OXi 

It is seen that the mass flux and the entropy flux have different velocities both in magnitude 
and direction. This means that a mass flux without entropy transport can take place. This 
implies an existence of a superfluid density p*. 
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In the case of quadratic dispersion law 

pii = mn^6ik, ^ = — , Pii = Su, ii = ViiKi (49) 
D n 

and the additional entropy flux vanishes. This implies that the superfluid effects should be 
negligible at very low temperatures where the excitations with parabolic dispersion relation 
predominate. 

Let us now rewrite the hydrodynamic set in terms of Landau superfluid theory in order 
to see better the analogy. 

Cubic crystal 

Let us flrst consider, for simplicity, a cubic crystal. Then, to the second order with respect 
to velocities one has 

J(r,t) =nV + z/W, Q(r,t) = iy°V+ (p + g)W (50) 

Iiiu=p5ik, L,, = -fi0 5,fc, F = 5V + j9tW (51) 



where 



2 

1 f I dek\ a „ m 



= (^ifj -^^k, ^ -^n^ (52) 

u{v,t) = ^rn j -gj^nldk, ^ —n (53) 
g(r,t) = im|..^<dk, ^ (54) 



and the following relations take place: 



'dp\ f dp\ p + q-fiu 



The meaning of the quantities involved can be seen from their limiting expressions in the 
effective mass (m*) approximation. 

In the notation of Landau theory: F = SV", V = V and the system of equations 
dM®]) takes the form: 

h + nMiv + n"div V" = 0, S + Sdiv V" = 0, (56) 

K + SVT + nV/i = 0, ptVT + p^V/i = (57) 



where — S—, — n — and the number density flux is 
Pt 



J = n'V + n"V" 



Second Sound 



One has in variables fj,, T, V*, 



aT + /3fi + n'V.y + n'^V.V" 
SVT + nV/i + p'Y' + p"V^ 



where 



9n 



dS_ 



P 



dn 
djjL 
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K = pV + mnW = p"V" + p"V" 
For quasiparticles with quadratic dispersion law — — 0: 



aT + /3/i + nV. V" 
,SVr + nVp + p'Y' + pV" 









If the number of quasiparticles is not conserved 



^o(q) = 



If the number of quasiparticles is conserved 



an 



TS 
Cvp 



(3S 










as; 



T / ds 



Cy \dvo 



dp 

dvo/ T 



p 



dP\ n 
dn J p^ 



12 



where s — S/n, Vq — 1/n . 



Quasiparticles with non-quadratic dispersion law 



If the number of quasiparticles is not conserved 



a;^(q) = (l + 5'')— -q^ = 




(67) 



If the number of quasiparticles is conserved 




(i-0^'(q,p= 



np" a {dP/dn)s ^ J 
')+5-u:l{c^,p = p-)-ujl{^) 



1 TS {dP/dT)l ^\ 



(68) 



P' 



n 



(69) 




0^2 (q) > cc;i(q) 



(70) 



Conclusion 

It is shown that the theory of superfluidity of solids should not be a replica of the Lan- 
dau theory of superfluidity. For crystalline bodies a two- velocity theory of supersolidity is 

presented with accounting of the quasimomentum conservation law. Such a theory cannot 
be applied to disordered systems, glasses etc. 
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